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We present a simple randomized scheme for triangulating a set P of n points in the plane, 
and construct a kinetic data structure "which maintains the triangulation as the points of P 
move continuously along piece'wise algebraic trajectories of constant description complexity. Our 

\^0 ' triangulation scheme experiences an expected number of (9(?7^/3s+2(^) log n) discrete changes, 

and handles them in a manner that satisfies all the standard requirements from a kinetic data 
structure: compactness, efficiency, locality and responsiveness. Here s is the maximum number 
of times "where any specific triple of points of P can become coUinear, /3s+2i<l) = ^s+2{<l) 1 1, and 

v»i«/ ' As+2('?) is the maximum length of Davcnport-Schinzel sequences of order s + 2 on ?7 symbols. 

c/3 , Thus, compared to the previous solution of Agar"wal ct al. [4], "we achieve a (slightly) improved 

bound on the number of discrete changes in the triangulation. In addition, "we believe that our 
scheme is simpler to implement and analyze. 
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1 Introduction 

Let P{t) = {pi{t), . . . , pn{t)} he a set of n moving points in the plane. We assume that the motions 
of the points are simple, in the sense that the trajectory of each point is a piecewise-algebraic curve 
of constant description complexity, meaning that it can be described as a Boolean combination of 
a constant number of polynomial equalities and inequalities of constant maximum degree. 

Our goal is to devise a reasonably simple scheme for triangulating P{t) at any fixed time t, and 
to maintain the triangulation as the points move. That is, we wish to partition the convex hull 
C'H{P) of P into pairwise openly disjoint triangles whose vertices are the points of P, so that the 
interior of each triangle is empty — it does not contain any point of P. The scheme has to be kinetic, 
so that we can keep track of the discrete combinatorial changes that the triangulation undergoes 
as the points move, and update the triangulation so that it continues to conform to the underlying 
scheme. (That is, at any given time t the maintained triangulation coincides with the one that 
would result in applying the static scheme to P{t).) 

The study of triangulations plays a central role in computational geometry because of their nu- 
merous applications in such areas as computer graphics, physical simulation, collision detection, and 
geographic information systems [8l[T3]. With the advancement in technology, many applications, 
for instance, video games, virtual reality, dynamic simulations, and robotics, call for maintaining 
a triangulation as the points move. For example, the arbitrary Eulerian-Lagrangian method [12] 
provides a way to integrate the motion of fluids and solids within a moving finite-element mesh. 

In M^, the Delaunay triangulation DT{P) of P produces well-shaped triangles, and thus is a 
good candidate for such a triangulation scheme. The problem, though, is that the best known 
upper bound on the number of discrete changes in DT[P[t)), as a function of time t, is only nearly 
cubic in n (the bound is cubic if the points move with constant velocities); see [2t ll4t [T5 l l20j. While 
it is strongly believed that the maximum possible number of discrete changes that DT{P{t)) can 
experience is only nearly quadratic in n, this is one of the hardest open problems in computational 
and combinatorial geometry (as recognized, e.g., in [H]). Until this conjecture is established, 
one seeks alternative triangulation schemes with a provable nearly- quadratic upper bound on the 
number of discrete changes. (This is best possible, since the convex hull itself can change O(n^) 
times during a simple motion of the points of P; see [20].) Moreover, the scheme should be 
sufficiently simple to define, to implement, and (as a secondary aesthetic virtue) to analyze. Finally, 
the scheme should satisfy the four basic properties of kinetic data structures [7] detailed below. 

Agarwal et al. [4] have recently presented such a randomized triangulation scheme which ex- 
periences 0(n^2^^^'^ ^^ s'^) discrete changes. Their scheme, however, is fairly complicated, and 
its analysis is also rather involved. It uses a hierarchy of subsets 9 = Rq (^ Ri Q • • • 'Z Rw = P, 
where each set Ri-i, for 1 <i <w,\s a. random sample of roughly |i?j|^~^'' logn points of Ri. The 
algorithm maintains an entire hierarchy of triangulations = 7o^7i^---^7^ = 7', where each 
7i is a triangulation of Ri; it is a refinement of 7i-i which is obtained by a suitable variant of the 
fan triangulation, introduced in [1]. 

Kinetic data structures. The Kinetic data structure (KDS) framework, introduced by Basch 
et al. [7], proposes an algorithmic approach, together with several quality criteria, for maintaining 
certain geometric configurations determined by a set of objects, each moving along a trajectory 
whose graph, as a function of time, is a piecewise-algebraic curve (in space-time) of constant 
description complexity. Several interesting algorithms have been designed, using this framework, 
over the past decade, including algorithms for maintaining the convex hull of a set of (moving) 



points in the plane [7] , the closest pair and all nearest neighbors in any dimension [3l [7] , and many 
other configurations. See [16] for a comprehensive, albeit old, survey, and [3] for a list of more 
recent results and references. 

Typically, a KDS operates by maintaining a set of certificates. As long as they are all valid, the 
structure being maintained is guaranteed to be valid too. Each certificate has a (first future) failure 
time, and we store these critical times in an event priority queue. When a certificate fails, we repair 
the KDS, update, if needed, the geometric structure that we maintain, generate new certificates 
and insert their failure times into the queue. 

Generally, a good KDS is expected to possess the following four properties: (i) Compactness, 
meaning that the storage that it requires is larger only by a polylogarithmic factor than the space 
required for the structure being maintained, (ii) Efficiency, meaning that the number of events 
that it processes (i.e., failure times of the certificates) is larger only by a polylogarithmic factor 
than the maximum possible number of discrete changes in the structure being maintained, (iii) 
Responsiveness, meaning that repairing the KDS at a certificate failure event takes only polylog- 
arithmic time, (iv) Locality, meaning that each input object is stored at only a polylogarithmic 
number of places in the KDS, so that an expected change in the motion of a single object can be 
processed efficiently. See O [7] for more details. 

Therefore, a good KDS for kinetic triangulation in M? should have only nearly linear storage, 
process only a nearly-quadratic number of events, each in polylogarithmic time, and each moving 
point should be stored at only a polylogarithmic number of places in the KDS. 

Our result. In Section [2l we present a simple triangulation scheme for a set P oin moving points 
in the plane. For the sake of efficient kinetization we make the scheme randomized, and assume a 
(natural) model in which the flight plans of the moving points are independent of the randomization 
used by the algorithm. The basic idea of the (static) triangulation is quite simple (some details are 
glossed over in this informal overview): We sort the points of P by their x-coordinates, split P at 
a (random) point p into a left portion Pl and a right portion Pr, compute recursively the upper 
convex hulls of Pl U {p} and of Pr U {p}, and merge them into the upper convex hull of the whole 
set P. 

This process results in a pseudo-triangulation of the portion of the convex hull of P lying above 
the x-monotone polygonal chain C{P) connecting the points of P in their x-order. Each pseudo- 
triangle is x-monotone, and consists of an upper base and of a left and right lower concave chains, 
meeting at its bottom apex. See Figure [1] for an illustration. A symmetric process is applied to 
the portion of the hull below C{P), by computing recursively lower convex hulls of the respective 
subsets of P. (In particular, we obtain a hierarchical representation oi C'H{P), similar to the one of 
Overmars and van Leeuwen [18] ; see also [5] . See [H [U [21] for additional applications of hierarchical 
pseudo-triangulations to kinetic problems.) 

To obtain a proper triangulation of (the convex hull of) P, we partition each pseudo-triangle 
T into triangles. We accomplish this in the following randomized incremental manner. We process 
the vertices of r (other than its apex and its leftmost and rightmost vertices) in order, according to 
the random ranks that they received during the first splitting phase, and draw from each processed 
vertex v a chord, within the current sub-pseudo-triangle r' of r containing v, which splits r' into 
two sub-pseudo-triangles. This process ends with a triangulation of r, and we apply it to each of 
the pseudo-triangles, to obtain the full triangulation of C'H{P). 

In Section [3l we prove that the expected number of events that can arise during the motion 
is 0{in? fis+2{n) log n) (with s and /? as defined in the abstract), and that the expected number of 



discrete (also called topological) changes caused in our triangulation by each such event is bounded 
by O(logn). 

In Section HI we show how to maintain this triangulation, as the points of P move, using a 
kinetic data structure that satisfies the criteria of [7], as listed above. There are several kinds of 
critical events we need to watch for, in which pairs of points are swapped in the x-order or triples of 
points become collinear. We process each event of the former type in 0(log n) expected time, and 
each event of the latter type in 0(log n) expected time, for a total of 0(n^ /3s+2(n) log^ n) (expected) 
processing time. Our implementation encodes the pseudo-triangulation as a treap on P |19j . 

The upper bounds that we obtain on the number of discrete events, and on their overall pro- 
cessing time, are slightly better than those of the scheme in [4j, and we believe that our scheme is 
simpler (and more "explicit") than that of [1]. 

2 The Static Triangulation 

In this section we describe a simple scheme for constructing a static triangulation T{P) oiCT-L{P). 
We fix a random permutation vr of the points of P. For each p G P we denote its rank in ir 
as priority (p). Let C{P) denote, as above, the x- monotone polygonal chain which connects the 
points of P in their x-order, assuming that no two points of P have the same x-coordinate. (In 
degenerate cases, which will arise at discrete instances during the motion of the points of P, C{P) 
connects the points in the lexicographical order of their coordinates.) Since the two points of P 
with extreme x-coordinates are vertices of CT-L{P), C{P) partitions C'H{P) into two components, 
C1-C^{P) and C'H~{P), lying respectively above and below C{P). With no loss of generality, we 
only describe a triangulation T~^{P) of C7i~^{P), and obtain the triangulation T~{P) of C'H~{P) 
in a fully symmetric fashion. The overall triangulation T{P) is the union of T^{P) and T~ {P). 

A static pseudo-triangulation o{CT-l'^{P). We first construct a pseudo-triangulation ofC'H'^{P) 
and then refine it into a triangulation by partitioning each pseudo-triangle into triangles. 

Each pseudo-triangle r that we construct consists of a left tail, a middle funnel, and a right 
tail (any of these substructures may be empty; the tails were not inentioned in the overview in 
the introduction). The funnel is an x- monotone simple polygon, whose boundary consists of an 
upper base, which is the segment connecting its leftmost and rightmost vertices, and of a left and 
right lower concave chains, which are denoted respectively as L(r) and 'R-{t). The point in which 
L(r) and 7^(t) meet is called the apex of r and denoted by apex(T). The left chain L(r) extends 
from the left endpoint of the base to apex(T), and the right chain extends from apex(r) to the 
right endpoint of the base; see Figure [TJ In addition, r may have a left tai^l L~ (t) and a right 
tail TZ~^{t), so that L"(t) is an x- monotone polygonal chain which extends from the left vertex of 
the funnel to the left, till the left endpoint left(T) of r, so that L~(t) U L(t) is a concave chain, 
and symmmetrically for TZ~^{t), which extends to the right till the right endpoint right(r) of r. 
Moreover, the line containing the base of r is an upper common tangent of L~(t) U L(r) and 
7^(t) L)TZ~^{t). Again, see Figured! We construct the pseudo-triangulation oi CH^iP) recursively. 
At each step of the recursion we have some subset Q C P of points which are consecutive in the 
x-order of P, and we construct a pseudo-triangulation VT^{Q) of C'H~^{Q). At the topmost level 
of the recursion we have Q = P. The construction of VT^{Q) proceeds as follows; see Figured! 
Let left((5) (resp., right(Q)) denote the point of Q with the minimal (resp., maximal) x-coordinate, 



^Tliese tailed pseudo-triangles are a special case of so-called geodesic triangles introduced in ^. 
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Figure 1: A single pseudo-triangle r in our pseudo-triangulation of CH^(P). In addition to its funnel (drawn 
shaded), t has two tails L^(t),7?.+ (t). 



and let m\d{Q) be the point p of Q\ {\eft{Q), nght{Q)} with the minimum value of priority(p). Set 
Ql = {p ^ Q \ x{p) < x{m\d{Q))}, Qr = {p € Q | x{p) > x(mid((5))} (so mid(Q) belongs to both 
sets). We add to VT^{Q) the following pseudo-triangle r. The base of r is the portion of the upper 
common tangent to C'H~^{Ql) and C'H^{Qr) between the points of tangency. We call this base the 
bridge of r and denote it by bridge(r). The left (resp., right) chain L(t) (resp., 7?-(r)) is the portion 
of the upper hull of Ql (resp., Qr) below bridge(r). We take L~ (r) to be the portion of the upper 
hull of Ql to the left of L(r), and define TZ'^{t) symmetrically as the portion of the upper hull of 
Qr to the right of 'R.{t). The points left((5) and right((5) become the respective endpoints left(r), 
right(r) of r. We also have apex(r) = rr\\d{Q) which belongs, by definition, to both chains. (The 
funnel of r may be empty, if mid(Q) is a vertex of the upper hull of Q. In this case we can think 
of the funnel of r as the singleton apex(r) = mid(Q), and r consists of the two tails L~{t),TZ~^{t), 
meeting at mid(Q), and forming together a common concave chain. Similarly, a pseudo-triangle 
may have an empty left tail and/or empty right tail.) 
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Figure 2: The recursive pseudo-triangulation of CH^{Q). We add to VT^iQ) the pseudo-triangle t (whose 
funnel is drawn shaded), with endpoints left(T) = left((5), right(r) = right(Q), and apex(r) = mid(Q), and then 
recursively construct VT^{Ql),'PT^{Qb)- 



We then recursively pseudo-triangulate each of CT-L~^{Ql), C'H~^{Qr). The recursion terminates 
when IQI < 3 (by construction, \Q\ > 2). If \Q\ = 3 then we output a single pseudo-triangle r, which 
is either a triangle, when the midpoint lies below the segment connecting the endpoints, or, in the 
opposite case, consists of the two segments L~(t) = left(T)apex(r) and TZ'^{t) = apex(r)right(T). 
If \Q\ = 2, no pseudo-triangle is output. In this case C7i^{Q) is a single edge of the chain C{P). 

Consider a pseudo-triangle r such that left(r) is not the leftmost point of P and right(r) is 
not the rightmost point of P. Then one can show that the triple (left(r), right(r), apex(r)) have 
the smallest priorities among all points whose ^-coordinates are between x(left(r)) and a;(right(r)), 
inclusive (see Lemma [2.11 below) . To make this true for all pseudo-triangles, we augment the initial 
point set P with two dummy points p-oo = {—oo, — oo) and Poo = (oo> —00)5 and assign to them 
priorities —1 and 0. The upper hull of the augmented point set is obtained from the upper hull of P 



by adding two vertical downward-directed rays at the leftmost and rightmost points of P. Hence, 
any triangulation of C7i^{P) is also a triangulation of C7i~^{P U {p^oo,Poo}), and vice versa. In 
the rest of the paper we denote by P the augmented point set. 

The following lemma gives an operational definition of VT^{P), which will be used in the 
sequel. 

Lemma 2.1. Let a, b, and c be three points in P, such that x{a) < x{b) < x(c). Then VT^{P) 

contains a pseudo-triangle r having endpoints left(r) = a, right(r) = c, and apex(r) = b, if and 

only if 

(i) priority(6) > max{priority(a), priority(c)}, and 

(a) all points p G P, such that x{a) < x{p) < x{c) satisfy priority (p) > priority (6). 

Proof. To prove the "only if" part we proceed by induction on our recursive construction. Recall 
that at each recursive step we process some subset Q ^ P whose points are consecutive in the 
x-order of P, and add to VT'^{P) a pseudo-triangle r with left(r) = left(Q), right(r) = right(Q), 
and apex(r) = mid(Q). To establish both asserted conditions (i) and (ii) for r, it is sufficient 
to observe that each point p, such that a;(left(Q)) < x{p) < a;(right(Q)), satisfies priority (p) > 
max{priority(left(<5)), priority(right(Q))}. Indeed, the desired property holds initially for P by our 
choice of the artificial points p-oo and poo and their priorities. Assuming that this holds when we 
process some subset Q, and using the fact that rr\\d{Q) is the point with smallest priority in the 
range a;(left(Q)) < x{p) < x{r\ght{Q)), the claim also holds for Ql and Qr. 

For the "if" part, we observe that for every choice of 6 G P there is exactly one choice of a and 
c in P so that the triple (a, b, c) satisfies (i) and (ii), and every point 6 € P is an apex of exactly 
one pseudo-triangle of 'PT^{P) (and the apex of each pseudo-triangle is distinct from each of poo 
and p-oo)- The latter is easy to establish by induction on the increasing order of the priorities of 
the points. This, combined with the arguments in the "only if" part, completes the proof. D 

The pseudo-triangulation tree. The pseudo-triangulation VT~^{P) can be represented by a 
binary tree in which every node v represents a pseudo-triangle Ty € VT^{P), and stores the 
point Pv = apex{Ty). The inorder of the tree is the increasing x-order of the apices (i.e., the 
points of P). The subtree rooted at v represents the recursive pseudo-triangulation of C'H^(Pu U 
{left(T^), right(T„)}), where Pv ^ P denotes the set of points stored at the nodes of the subtree 
rooted at v. Note that \eft{Ty) and right(T^,) are not stored at this subtree — they are the next points 
to the left and to the right of the points of P„. Abusing the notation slightly, we denote by VT~^{P) 
both the pseudo-triangulation VT'^iP) and the tree representing it. 

Remark: Let v be a node in VT'^iP), so that left(r^) ^ p^oo- Then left(rt,) is stored at the lowest 
ancestor of v whose right subree contains v. If left(T^) = p-oo then v belongs to the path from the 
root of T'T~^{P) to the leftmost leaf. Symmetric properties hold for right(T^). 

In summary, we have the following lemma, whose proof is immediate from the construction. 

Lemma 2.2. The tree representing VT^{P) is a treap on P\ {p-ooPoo}- That is, VT^{P) is a 
heap with respect to the priorities, and a search tree with respect to the x-coordinates of the points. 

Triangulating a fixed pseudo-triangle. Let r be a pseudo-triangle of VT'^{P). Assume that 
the funnel of r is not empty, and is not already a triangle. We say that two vertices p, q of the 
funnel of r, where p belongs to L(r) and q belongs to 7^(t), are visible from each other if pq does 
not intersect dr (except at its endpoints); in this case pq lies inside the funnel of r. Denote by 



^{p) the rightmost point on the right chain which is visible from p. Note that either i'{p) is the 
rightmost vertex of r or py{p) is an upper tangent to TZ{t). Symmetic definition and properties 
hold for points q on TZ{t). This definition also applies when p is the leftmost vertex of L(r) and 
when q the rightmost vertex of TZir) (the endpoints of bridge(r)), in which case z/(p) = q and 
^(q) = P- See Figure [3] (left). 
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Figure 3: Left: The first step of triangulating a single pseudo-triangle r S 7^7~~''(P). Right: During the recursive 
construction of T{t) every non-corner vertex p of the funnel of r generates exactly one edge Cp = pv*{p), thus 
recursively splitting some sub-pseudo-triangle Tp (drawn shaded). Note that in this figure v*{p) ^ v{p), which is 
the left endpoint of bridge(T). 



The triangulation T{t) of r is obtained by recursively splitting r by chords into sub-pseudo- 
triangles, in the following manner. Choose the minimum priority vertex q of the funnel of r, other 
than the leftmost and the rightmost vertices and the apex. Assume, without loss of generality, 
that q lies on L(r). See Figure [3] (left). The segment qui^q) splits r into two sub-pseudo-triangles 
r"*" and r~. The pseudo-triangle r"^ has q as an apex and the same base as r. Its left chain is 
the portion of L(r) from q to the left, and its right chain is the concatenation of qy{q) with the 
portion of TZ{t) to the right of i'{q). The pseudo-triangle r^ has qi'{q) as its base, the same apex 
as r, and its left and right chains are the portions of L(r) and TZ{t) delimited respectively by q 
and by ^{q). A symmetric situation arises when q G "^(t). We add the edge qv{q) to T{t), and 
recursively triangulate each of r"*" and r~. We say that the edge qi'{q) in T{t) is generated by q. 
In the further recursive steps, we redefine v{p), for vertices p of each of these sub-pseudo-triangles, 
restricting the visibility to only within the respective pseudo-triangle. Note that for any pair of 
vertices p, q that lie on the same chain of r, the segments pi'{p) and qi'{q) do not intersect in their 
relative interiors. Therefore, if ^{p) changes after a recursive call then it must change to a vertex 
of the base of the corresponding sub-pseudo-triangle. See Figured (right). The recursion bottoms 
out when the interior of r is a triangle. Note also that all the chords in T{t) cross the vertical ray 
above apex(r), and so they are totally ordered in the vertical direction. 

Properties of T{t). Every vertex p of the funnel of r, other than the leftmost and the rightmost 
vertices and the apex, generates exactly one edge Cp during the whole recursive process. (For 
example, in Figure [3] (left), the vertex i/(g) will not generate an edge in r~, since it is an endpoint 
of that funnel, but will still generate an edge within r"*", or within some recursive sub-pseudo- 
triangle of r"*".) We denote by Tp the sub-pseudo-triangle in which e^ is generated, and by i'*{p) 
the other endpoint of Cp. Note that i'*{p) is either the original z/(p) or an endpoint of the base of 

T„. 



3 Number of Discrete Changes in T{P) 

In this section we bound the overall expected number of discrete changes that T{P{t)) experiences 
as the points of P move along (continuous) pseudo-algebraic trajectories of constant description 
complexity. The analysis is with respect to a fixed random permutation vr of P drawn ahead of the 
motion, so that the motion is "oblivious" to the choice of vr. Thus, even though the x-order of the 
points may change during the motion, each point retains its initial priority, and the permutation vr 
is still a random permutation of P, with respect to the x-order of these points, at any fixed t. 

Discrete changes in 'PT^{P). For a fixed time instance t E M, each pseudo-triangle r G 
'PT'^[P{t)) is defined by its endpoints left(r), right(r), and by its apex apex(r). Given such a 
triple of points, they define a valid pseudo-triangle at time t if and only if they, and the points 
in-between in the x-order, satisfy the conditions of Lemma |2. II (at time t). Thus, as long as the 
x-order of the points does not change, VT^{P{t)) does not change either. That is, it consists of a 
fixed set of pseudo-triangles, each defined by a fixed triple of points. However, the geometric struc- 
ture of a pseudo-triangle may change during such a time interval, and we will bound the number 
of these changes separately. Changes in (the labelings of the pseudo-triangles of) VT~^{P{t)) occur 
only at discrete times when the x-order of some pair of points in P{t) changes; we refer to these 
changes as x-swap events. 

We assume that each pseudo-triangle r is present in VT'^{P{t)) at a maximal connected time 
interval /(r), which is associated with r. That is, pseudo-triangles with the same triple left(r), 
right(r), and apex(T) that appear in VT~^{P{t)) at disjoint time intervals, are considered distinct. 
We emphasize that all the other features of r, such as bridge(r), the chains L(r) and TZ{t), and the 
triangulation 7~(t) of its funnel, may undergo discrete changes during the time interval /(t). A 
pseudo-triangle r is created or destroyed only at a swap event when a point p £ P with priority (p) < 
priority(apex(r)) crosses one of the vertical lines through its endpoints left(r) and right(r) (of course, 
this also subsumes the cases where priority (p) is smaller than that of an endpoint of r), or when 
the x-order of the points in the triple defining r changes. In the former case, if priority (p) > 
max{priority(left(r)), priority(right(r))} then r is replaced by another pseudo-triangle r' with the 
same endpoints left(r') = left(r), right(r') = right(r) but with p as a new apex. 

If priority(p) < max{priority(left(r)), priority (right(r))} then p replaces the endpoint it was swapped 
with. Thus, each pseudo-triangle r in our kinetic pseudo-triangulation 

VT~^{P{t)) is defined by at most five points: left(r), right(T), apex(r), and at most two additional 
points which determine, by swaps with the endpoints of r, the endpoints of the lifespan I{t) of r 
inPr+(P(t)). 

Discrete changes in T{t). Fix a pseudo-triangle r G VT^{P{t)). We only consider discrete 
changes in the funnel of r and its triangulation /(t), and ignore changes in the tails L~(r), 7?.~''(r) 
(unless they also affect the funnel). This is because the changes in the tails will also show up as 
changes in the funnels of other pseudo-triangles that are created further down the recursion. 

For a fixed time instance t € I{t), the combinatorial structure of the triangulation T{t) of r 
depends only on the discrete structure of the boundary of the funnel of r (i.e., the ordered sequences 
of the points along the chains L(r), 7^(t), and the base bridge(r)) and the visibility points v^p) of 
all the vertices of the funnel of r, excluding apex(r) (of course, it also depends on vr). Therefore, 
as the points of P move during the time interval /(r), 7~(t) can change combinatorially only at 
events where the boundary or visibility structure of r changes. These events fall into the following 
three types: 

7 



(i) Envelope events, which occur at instances when one of the chains L(t), 7^(t) contains three 
colhnear vertices; see Figure H] (right). This happens when a vertex (which is not an endpoint of 
bridge(r)) is added to or removed from one of the chains bounding r. We denote the total number 
of such events during the period /(r) by Er- 
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Figure 4: Left and right: visibility and envelope events (respectively). The sub-pseudo-triangle tq contains all 
edges which are inserted to or deleted from T{t) at this event. 

(ii) Visibility events, at which a vertex q of TZ{t) becomes coHinear with an edge pr of L(r) U 
L^(t), or vice versa. See Figured] (left) (L~(t) is relevant only for visibility events that affect 
bridge(T) and then pr has to be its rightmost edge, and symmetrically for TZ'^{t)). This happens 
when ^{q) changes from p to r, or vice versa. In particular, each (discrete) change of bridge(r) 
corresponds to a visibility event in which the bridge becomes collinear with an edge of L~(t) UL(r) 
or of TZ{t) U TZ'^{t) that is incident to the respective endpoint of the bridge. We denote the total 
number of visibility events during /(r) by Vr- 

A special case of this event occurs when bridge(r) is created (resp., destroyed), so that right 
before (resp., after) the event, the funnel of r is empty. Note that immediately after (resp., before) 
the creation (resp., destruction) of bridge(T), the funnel of r is a triangle. 

(iii) Swap events, at which some point p G P, satisfying priority (p) > priority (a pex(r)), crosses 
one of the vertical lines through left(r), right(T) or apex(T). Note that a single swap event of this 
kind may cause massive discrete changes, of highly unlocal nature, in the chains L(r), TZ{t), in the 
visibility pointers v{q) of the vertices of r, and in bridge(T). See Figured] (left) for an illustration. 

Note that a swap between any other pair of points p, q within the x-range of r can be ignored 
in the present analysis, since the lower of the two points cannot belong to the funnel of r at the 
time of swap. 

Assuming general position of the trajectories of the points, the above events occur at distinct 
time instances (except that the same event may show up, in different forms, in several pseudo- 
triangles). 

A visibility event happens when v{p) changes for some point p; we then say that p is involved 
in the visibility event. An envelope event happens when a point p joins or leaves one of the chains 
L(r), TZ{t); we then say that p is involved in the envelope event. 

Lemma 3.1. The only point p for which ^{p) changes in an envelope event is the point p involved 
in the event. 

Proof. The lemma follows since at the moment following (resp., preceding) the appearance of p on 
(resp., disappearance from) its chain, say L(r), its two incident edges are almost collinear. Thus, 
all vertices q on the opposite chain satisfy u[q) ^ p both before and after the event, and u[q) is not 
affected by the event. D 



The number of changes in T{t). We define Pr as the set of points p ^ P, other than apex(r), 
that appear on C{P) between left(r) and right(r), at any time during the life span /(r), and put 
Nr = \Pt\- (Note that the points of Pr may enter or leave the interval between left(r) and right(r) 
in the middle of /(r), at x-swaps with either left(r) or right(r).) As noted above, every point 
p & Pr satisfies priority(p) > priority(apex(r)). Clearly, our triangulation undergoes 0{Nr) swap 
events during /(r) (recall that we only consider swaps with left(r), right(r) or apex(r)), and each 
of them leads to 0{Nr) edge insertions and deletions to T{t) (the maximum number of edges in 
the whole triangulation T{t)), for a total of 0{Nr) such updates. We next bound the number of 
discrete changes in T{t) caused by events of the remaining two types. 

Fix a set of at most five points that can potentially define a pseudo-triangle for some set of 
priorities. This set has an associated time interval [ti,t2]; and consists of three points o, 6, and c, 
such that, at all times ii < t < t2-, x{a(t)) < x{b{t)) < x{c{t)), and of two additional points di and 
d2 (each of which could be equal to 6), so that the x-coordinate of di swaps with either a or c at 
times ti, for i = 1,2. For some drawings of the random priorities, r appears as a pseudo-triangle, 
and for other drawings it does not. For r to appear in VT^{P{t)), the priorities of a = left(r) and 
c = right(r) should be smaller than the priority of b = apex(r). The priorities of di and ^2 have to 
be at most the priority of 6 = apex(T), and the priorities of all other points in Pr should be larger 
than the priority of b. The probability of this to happen, assuming a, b, c, di, d2 are all distinct, is 
easily seen to be 0{l/N^) (for Nr > 0). 

When we condition on drawings in which r indeed appears in VT^{P(t)), the following holds. 

Proposition 3.2. Let t be a pseudo-triangle in the kinetic triangulation VT~^{P{t)). Then the 
expected number of discrete changes in the triangulation T{t) of r, after any single envelope or 
visibility event which happens during the period I{t), and conditioned on t appearing in VT^{P{t)), 
is 0{logNr) = 0{logn). 

Proof. Clearly, the chords of T{t) (the additional edges which partition r into triangles) admit 
a total vertical order, because they all cross the vertical line through apex(r). Consider a time 
instance to ^ I^ when an envelope or a visibility event occurs, and let po € L(r) UTZ{t) be the point 
involved in the event. Let t^ (resp., t^) be the time right before (resp., after) the event. Note that 
Pq cannot be the apex of r (unless the funnel of r is already, or is going to become, a triangle). Note 
also that po is not a vertex of bridge(r), neither at t^ nor at t^, unless po is involved in a visibility 
event which changes bridge(r). In the latter case, T{t) gains or loses its topmost triangle at time to 
but there are no other changes in the triangulation, as is easily checked. We may therefore assume 
that bridge(r) does not change at time to, and that po is not a vertex of bridge(T). 

With no loss of generality, we assume that pQ is a vertex of TZ{t) at time t^, and treat the 
remaining cases symmetrically (for a visibility event, po belongs to TZ{t) also at time t^). Consider 
the triangulation T{t) at time t^ (that we would have obtained if we were to reconstruct T^{P) 
statically at time t^). Let tq be the sub-pseudo-triangle of r within which the edge Po^*(Po) is 
generated during the construction of T'^{P) (see Figure H] (right)). Note that the event at time to 
leaves unchanged the visibility vertex u[p) of each vertex p in r other than pQ. Indeed, this follows 
from Lemma l3. II for envelope events and is obvious for visibility events, using our assumption that 
bridge(r) does not change. The recursive construction of T{t) is easily seen to imply that tq appears 
as a sub-pseudo-triangle in the construction also at time t^. Indeed, an easy inductive argument 
on the order of the ranks of the funnel vertices implies that the modified visibility vertices v* {p) , 
and the resulting chords pv*{p), also do not change, up to the point where tq is constructed. Right 



after this step, the chord from po is drawn, so the rest of the construction of Tir) might change 
completely, but only within tq. Hence, tq contains every edge which is inserted to or deleted from 
T{t) at time to. Therefore, the number of changes in T{t) is bounded by 0{Wo), where Wq denotes 
the number of vertices of tq at the time of the event. 

Note that Wq is a random variable depending (only) on the permutation 7r{Pr) of P^-, which 
is obtained by restricting vr to Pr- Recall that we condition the analysis on permutations vr such 
that r indeed appears in VT^{P). In these permutations, the points of Pr have to follow all the 
(at most) five points defining r, but as long as they obey this restriction they can appear in any 
order. It follows that, in our conditional probability subspace, the restriction of vr to P^ is a random 
permutation of P-r- 

To bound the expected value of VFq, we fix an arbitrary threshold /c > 10 and prove that the 
event {Wq > k} occurs with probability at most 0{l/k). The expected value of Wq is then bounded 
by 

log Nr /log Nr \ 

J2 2'^'Pr{Wo>T} = oi Y, 1 =0(logiV.)- (1) 

i=0 V i=0 J 

To show that Pr{Wo > k} = 0{l/k), we proceed through the following cases. In each case, 
except for the last one, we find a set 5o of ^{k) points which does not depend on '7T{Pr), so that 
all its elements must appear in 7r{Pr) after pQ. This readily implies the asserted bound. The last 
case is more involved but it is still based on the same general idea. 

Visibility event. If v*{po) is a vertex of the base of tq, both at time t^ and at time t^, then T{t) 
does not change combinatorially at time to- Otherwise, as follows from the discussion in Section O 
all three vertices that become collinear in the event appear in tq, both before and after the event, 
which implies that i'{po) = i^*{po) at both times tQ and t^ (although they assume different values 
of these times). 

Recall that p is assumed to be a vertex of 7?.(r), and suppose that Wq > k. If tq contains at 
least k/2 vertices of '/^(t), then it also contains a sequence 5o of k/A — 1 consecutive vertices of 
TZ{t) either immediately to the left or immediately to the right of po- Otherwise, tq contains i^(po) 
together with at least k/2 — 1 other vertices of L(r), so it must contain a sequence Sq of k/4 — 1 
consecutive vertices of L(r) lying either immediately to the left or immediately to the right of i'{pq). 
In both cases, the key observation is that Sq does not depend on vr(PT-), and that po precedes all 
the vertices of 5*0 in 7r{Pr) (except possibly for one extremal vertex which is a corner of tq). As 
noted above, this establishes the asserted bound. 

Envelope event. Again, suppose that Wq > k. If tq contains at least k/2 vertices of 7^(t), the 
bound follows by exactly the same argument as in the case of a visibility event. Otherwise, if tq 
contains apex(r) we set Sq to be the first k/2 — 2 points of L(r) to the left of apex(r). Again, 5*0 
does not depend on 7r(PT-), and all its elements must appear in ii{Pt) after pQ, so the bound follows. 
We therefore assume that tq contains at most k/2 vertices of 7^(t), and that its apex go is 
distinct from apex(r). Thus, the edge qoiy*{qo) that go generates is the lowest edge of tq which is 
a chord of r. We argue that i'*{qo) = i^iQo) (before and after to; the definition of tq implies that 
qo precedes po in 7r(P,-)). Indeed, otherwise, by the definition of T{t), u*{qo) is a vertex of the 
base of tq, which happens only if one of the chains of tq consists of the single edge qov*{qo). Since 
Po £ T^ij) and is involved in an envelope event, the edge qov*{qo) must be the only edge of the 
left chain of ro, which contradicts the fact that L(t) must contain at least k/2 vertices of To (for 
k > 10). We distinguish between the following two cases. 
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(i) qo € L(t) (as depicted in Figure H] (right)). Then the entire left chain of tq is contained in 
L(r). Let io be the Hne passing through pQ and the other two vertices of 7^(r) participating in 
the envelope event, and let cq be the edge of L(r) intersected by £. Clearly, cq is contained in tq, 
because otherwise TZ{tq) would not be convex. If tq contains k/4 — 1 consecutive vertices of L(r) 
which lie immediately to the left cq, we set Sq to be the set of these points, except for the leftmost 
one (which may be the endpoint of the base of tq). Otherwise we set So to be the set of k/A — 2 
points lying on L(t) to the right of eg. Since the definition of cq does not depend on 7r(PT-), the set 
So too does not depend on n{Pr)- 

(ii) qo G ^(t) (as depicted in Figure [5] (right)). In this case we define at most k/2 sets, each 
consisting of n{l/k) points and independent of tt{Pt-), such that all the points in at least one of 
these sets appear after both po and qo in ■k{Pt)- We fix qo on 7^(r) to the left of po and define Sq^ as 
the set of k/2 — 2 consecutive vertices of L(r) which appear at time to (along L(t)) immediately to 
the left of i'*{qo) = ^((^o)- By the current assumptions, if qo is indeed the apex of tq then all points 
q S SqQ belong to tq and, hence, satisfy priority(q') > priority(po) > pi'ioi'ity(go)- Since qo is fixed, 
iSgQ is also fixed and is independent of i^{Pt)- Hence, the above event happens with probability 
0(l//c^). Moreover, go is one of the at most k/2 vertices of 7^(t) that lie to the left of po- Hence, by 
the probability union bound, the total probability of this scenario (over all the appropriate vertices 
qo G 7^(r)) is 0{l/k). 





Figure 5: Left: Swap event. The funnel of r immediately before the x-swap between p and apex(T), which 
causes the vertices pi and p2 to appear on L(t), and the vertices (71,(72, 93 to disappear from TZ{t). Right: 
Envelope event. The case in which go lies on TZ{t). 



We have proved that Pr (Wq > k) 
completes the proof of Proposition 



0{l/k) for any k > 10. This implies Equation ([T]) and 

D 

Corollary 3.3. Let t be a pseudo-triangle in the kinetic pseudo-triangulation VT~^{P{t)). Then 
the expected number of edge insertions and deletions in T{t) during the period I{t), conditioned 
upon the event that r appears in VT~^{P), is 0{{Et- + Vr) log A'^,- + N^). 

For a fixed pseudo-triangle r (including the choice of the connected life span I{t)), Vr and 
Ej- are 2-valued random variables: They are if r does not appear in VT~^{P), and assume a 
fixed "deterministic" value if r does appear. The following theorem gives an upper bound on these 
values. 

Theorem 3.4. For each pseudo-triangle r we have Vr = 0{N^I3s+2{Nr)) and Er = 0{N^I3s+2{Nr)), 
where s is the maximum number of times at which any fixed triple of points of P becomes collinear. 

Proof. We show the bound for visibility events. The bound for envelope events is known (see [HIS]) 
and can be proved similarly. 
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We fix a point p ^ P^- and count the number of visibility events where p is a vertex of L(r) 
which is coUinear with an edge of 7^(t). To do so, we define, for each q € P!^ = Pt-L) {right(r)} \ {p}, 
a partiahy defined function (pp^q{t) which measures the angle between pq and the y-axis, and whose 
domain consists of all t € M at which 3;(left(r)) < x{p) < rc(apex(T)) < x{q) < 3;(right(T)). Clearly, 
each visibility event under consideration corresponds to a breakpoint of the lower envelope of 
{'^P,q}qeP' (but not necessarily vice versa; for example, such a breakpoint can arise at a time when 
p is not a vertex of the funnel of r). Since any pair v^p^g^, V'p,g2 o^ these functions can intersect in 
at most s points (these are times at which p, qi, and q2 are collinear), and for each q the domain 
of ipp^q{t) consists of a constant number of intervals (delimited by times at which either p or q swap 
with left(T), right(r), or apex(T)), it follows that the number of breakpoints is 0(A'r/3s+2(^r)) J20] . 
A symmetric argument holds for the number of visibility events where p is a vertex of TZ{t) which 
is collinear with an edge of L(r). Repeating this analysis for each p ^ P yields the asserted overall 
bound. D 

Fix a pseudo-triangle r. Conditioned on priorities that cause r to appear in VT~^{P{t)), 
Corollary 13.31 and Theorem 13.41 imply that the expected number of discrete changes in 7~(t) is 
0{N^(3s+2{N^)logN^). Let Pr(T) be the probability that r indeed appears in VT^iPit)). Then 
the total expected number of discrete changes in VT^{P{t)) is 

O I Y.PriT)N^I3s+2{Nr) log Nr) =ohs+2{n)lognY,Prir)NA . 

Lemma 3.5. X^,- Pi'('7")-^r = O(n^logn), where the sum is over all (possible sets of 1 < h < 5 
points defining) possible pseudo-triangles r. 

Proof. Without loss of generality, we only consider pseudo-triangles r with A',- > 0, which are 
defined by five distinct points of P \ {p-ooPoo}- (Pseudo-triangles defined by fewer than five 
distinct points, or those whose defining 5-tuple includes p-oo and/or poo are analyzed similarly, 
replacing the exponent 5 by the appropriate 1 < /i < 4.) Thus, as already noted, Pr(r) = 0{1/N^), 
because r appears in VT'^{P{t)) if and only if the priorities of the five points that define r are 
smaller than the priorities of all other points in P^ (and apex(T) has the largest priority among the 

defining points). Therefore ^ Pr(T)iV^ = O ( ^ l/N^ | . 

In what follows, we call N^- the level of r. Let Mk{n) (resp. M<kin)) denote the maximum 
number of pseudo-triangles of level k (resp., of level at most k), defined by 5 points, in a set of n 
moving points. We claim that M<k{n) = 0{n'^k^). To see this, consider all the pseudo-triangles r 
(defined by five points) whose birth time is determined by a fixed x-swap event occurring at some 
time io; between some pair of points a, 5 € P. Assume without loss of generality that a = left(T). 
Then apex(r) and right(T) are among the k+2 points whose x-coordinates lie at time to immediately 
to the right of x{a) = x[h). Similarly, the fifth point, which is responsible for the destruction of 
r, is one of the first k + 1 points whose x-coordinates enter the interval between x(left(T)) = x(a) 
and x(right(r)). Thus, each of the 0{n'^) x-swap events defines the creation time of at most 0{k^) 
pseudo-triangles of level at most k, which readily implies the asserted bound on Ai<k{n). We thus 
have 

5^Pr(r)A2 = O [Y^l/Ni] = O [Y,M,{n)/e 

T \ r / \k>l 
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\k>l J \k>l J 

D 

The combination of Corollary 13.31 Theorem 13.41 and Lemma [3.51 implies the following summary 
theorem. 

Theorem 3.6. The total expected number of discrete changes in the kinetic triangulation T{P{t)) 
is O {n? I3s+2{n) \og^ n) . 

4 Kinetic Maintenance of T{P) 

In this section we describe a kinetic data structure which supports efficient maintenance of T^(i-'(t)) 
under motion. The structure satisfiecl the standard requirements of efficiency, compactness, re- 
sponsiveness, and locality, as reviewed in the introduction. 

The static structure. We store the pseudo-triangulation tree VT^{P) as a treap over P, as 
described in Theorem 12.21 whose inorder is the x-order of the points and where the heap order is 
according to their random priorities. Each node v in VT~^{P) corresponds to the pseudo-triangle 
Ty whose apex is the point stored at v. We also store at v, as auxiliary data, the endpoints left(rt,) 
and right(rt,), which are inherited from appropriate ancestors of v. 

In addition, we also store at v the combinatorial description of the funnel of r^,, and of its 
triangulation T{Ty). This includes bridge(r^), two ordered lists storing the vertices of L(t„), and 
TZ{Ty) in their left-to-right order, and the list of the chords of T{Ty), sorted in their vertical order 
(i.e., the order of their intersections with the vertical line through apex(T„)). We represent any 
sorted list of vertices or edgecl as a balanced binary tree supporting each of the operations search, 
split, and concatenate, in O(logn) time |22j . To facilitate efficient kinetic maintenance of T'''(-P(t)), 
we also store the vertices of the upper hull of P, in their left-to-right order in a balanced search 
tree. Note that each edge of the triangulation (not on C{P)) appears twice in our structure, once 
as bridge(Ti,) for some pseudo triangle r„, and once on L(r„,) or TZ{tw) for some ancestor w of v or 
on the convex hull of P. 

Theorem 4.1. Let P be a set of n points in the plane. The pseudo-triangulation tree VT^{P), 
augmented with the auxiliary data items, as above, uses 0{n) space, and it can be initialized in 
O(nlogn) time. 

Proof. The asserted bound on the overall storage follows from the easy observation that VT~^{P) 
contains 0{n) nodes, and every point p £ P appears as a non-corner vertex on at most one chain 
L(r^), TZ{Ty), over all nodes v of VT'^{P). 

We construct the pseudo-triangulation tree VT^{P) (excluding the auxiliary items bridge(rt,), 
L(r^), 7^(t„) and the chords of T{Ty)) in a single top-down pass, which implements the recursive 
construction given in Section [2l Clearly, this can be done in 0{n) time, after an initial sorting of 



^As in [5], all properties (except for compactness) hold in expectation, with respect to the random per- 
mutation TT. 

^Note that we do not store explicitly the tails L~(t),7^+(t), because the overall storage that they would 
require could be too large, as they can be shared by many pseudo-triangles. 
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the points of P, by their x-coordinates and by their priorities; sorting the points takes 0{nlogn) 
time. 

We next compute the items L(r„), TZ{Ty), and bridge(r^) stored at the nodes v of VT^{P), by 
a single bottom-up traversal of 'PT~^{P), which computes for every node v the upper hull Iy({v) of 
the set P^ U {left(rt,), right(ru)}. When we process a new non-leaf node v, we have already visited 
its respective left and right children V£ and Vr, so their hulls U{ve) and U{vr) are already available. 
We compute bridge(r„) in O(logn) time by a simultaneous binary seach over hl{vi) and V({vr), in 
the manner described in [H]. Then we use bridge(T^) to split ZY(f^) (resp., U{vr)) into L~(r„) and 
L(t^) (resp., TZ{Ty) and TZ~^{Ty)). We store explicitly the chains L{Ty) TZ{t^) at v, and compute U{v) 
by concatenating the three edge lists L~(r), {bridge(r)}, and 7^^(r), in a similar manner to that 
described in [5]. Overall, we spend O(logn) time at each node of 'PT'*' (-P) , for a total of O(nlogn) 
time. 

Finally, for each node v in VT'^{P), we compute the list of chords of T{Ty) using the recursive 
mechanism described in Section [2l Recall that every non-corner vertex p of the funnel of Ty 
generates exactly one edge Cp which recursively splits the unique sub-pseudo-triangle Tp of r^, . We 
process the non-corner vertices of T{Ty) in the increasing order of their priorities, and store the 
edges constructed so far in a list, in the order of their intersections with the vertical line through 
apex(T„). 

It takes O(logn) time to process a non-corner vertex p of r^,, for a total of O(nlogn) time. 
Indeed, we can determine the corners of Tp in O(logn) time, by a binary search over the list of 
the previously generated edges. In addition, we can determine i'(p) by a binary search over the 
appropriate chain L(r.i,) or TZ{Ty), obtain iy*{p) in 0(1) additional time, and insert the chord pv*{p) 
into the list of chords in O(logn) time. D 

The kinetic certificates. To ensure the validity of VT^{P) and its triangulation T'^(P), we use 
three types of certificates, denoted as CT, CE and CV. Each certificate is a predicate on a constant 
number of points. As long as all the certificates remain true, the validity of VT^{P) and T~^{P) 
is ensured. Each certificate contributes a critical event to the global event priority queue Q, which 
is the first future time at which the certificate becomes invalid (if there is such a time). 

CT -certificates. To ensure the validity of the tree VT^{P) (ignoring the auxiliary data), each 
pair of points p,q G P with consecutive x-coordinates contributes a CT-certificate asserting that 
the order of x{p) and x{q) remains unchanged. This certificate fails at the first future moment of 
an x-swap between p and q. According to Lemma 12.11 CT-certificates (together with the chosen 
priorities) are sufficient to ensure the validity of the "bare" tree VT'^{P). 

C^- certificates. For each node v in VT'^iP), the edge bridge(Tt,) = pq contributes a CE-certificate 
ensuring that the (current) neighbors of p and q on L^fr^,) U L(r^) and TZir^) U 7^+(r^,) remain 
below the line through p and q. This certificate involveqj at most six points and fails at the first 
future time of collinearity between p, g, and one of their four neighbor vertices on L~(rt,) U L(rt,) 
and on 7^(t^) U lZ^{Ty). 

So far, we have ensured the validity of the tree VT^{P) and of the edges bridge(rt,) stored at 
its nodes v. Moreover, the validity of all the chains L(t„),7^(tj,) is also ensured because each one 
of their edges either belongs to C{P) or appears as bridge(ri„) at some descendant w of v. Here 
a collinearity between three consecutive points on L(r„) or on TZ{Ty) (an envelope event) will be 



*If bridge(Tu) does not exist then we have an even simpler certificate which fails when the two edges of 
L~(tu),7?.+ (tu) incident to apex(Tu) become collinear. 
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detected as a change in bridge(r^), for the appropriate descendant w. Similarly, the validity of the 
upper hull of P follows since each of its edges either belongs to C{P) or appears as bridge(rt,) at 
some node v. See [5] and [18] for more details. 

C\/ -certificates. It only remains to ensure the validity of the triangulations T{Ty), over all nodes 
V € VT^{P)- For this we need the third type of certificates, denoted by CV. Fix a node v in 
VT^{P). Every internal point p of L(t^,) or TZ{ti,) contributes a CV certificate ensuring the validity 
of I'ip). This certificate involves p, v{p)-, and the two points adjacent to y{p) on its chain. It fails 
when one of the points adjacent to z^(p) becomes collinear with p and ^{p). 

Clearly, all of the above certificates use 0{n) storage, and can be initialized, including the 
construction of the event queue Q of their first failure times, by the algorithm of Theorem 14. !( 
without increasing its overall assymptotic running time, i.e., in 0(n log n) time. 

Handing critical events. We next describe the repair operations required when an event, at 
which some certificate fails, happens. 

QT -certificates. Failure of a CT-certificate occurs at an x-swap. That is, the order of the x- 
coordinates of two consecutive points along C{P) switches, at some time t = to- 

With no loss of generality we assume that priority(p) < priority(g'), implying that (/(tg ) is a 
descendant of p(tQ), where t^ ,1^ denote the time just before and just after t^, respectively. To 
update 'PT'^{P{tQ)) we reconstruct from scratch the subtree rooted at the node v containing p, 
and recompute the kinetic certificates associated with its nodes and the points that they contain. 
We remove the failure times of the expired certificates from Q, and insert the new ones. All this 
can be done in 0(n^logn„) time using the algorithm of Theorem 14.11 where n^ = \Pv\- We prove 
that E{n„} = 0(log?i) by applying a simplified version of the analysis used in Proposition 13.21 As 
above, it suffices to show that Pr{n^ > k} < A/k, for any A; > 1. Indeed, Uy > k implies that 
either each of the k/2 points w whose x-coordinates immediately precede x{p) or each of the k/2 
points w whose x-coordinates immediately follow x{p) at time to satisfies priority(p) < priority (tt;). 
This happens with probability at most A/k. Thus, we can reconstruct the subtree rooted at v in 
O(lognloglogn) expected time. 

As can easily be checked, if neither p nor q is the leftmost or the rightmost point of P (excluding 
the points at infinity which we added) then no further updates outside the subtree of v are needed, 
and no additional certificates need to be created or destroyed. (That is because left(ri,) 7^ P-00 and 
right(T„) 7^ Poo, so the upper hull U{v) contains at most one of p,q, and it does not change as a 
result of the swap.) We next describe the necessary modifications in the setting, depicted in Figure 
[H in which case we assume that (i) p and q are the two points with the smallest x-coordinates, (ii) 
x(g(tQ)) > x(p(t^)), and (iii) the y-coordinate of 7? is larger than at q when they swap; the other 
cases are treated symmetrically. The x-swap between p and q causes q to appear on the upper hull 
of P, below and to the left of p. We add q to the upper hull in O(logn) time. Similarly, q becomes 
part of the tail L~(t^) of every ancestor w of v (both w and v lie on the leftmost path of the treap). 
If w is such an ancestor whose bridge is incident to p (from the right), then we have to incorporate 
q into the certificate of bridge(T^), and possibly replace its old failure time in Q with a new one. 
Since the expected number of ancestors w of v, in the treap VT^{P), is O(logn) (see, e.g., [IH]), 
any swap event can be processed in 0(log n) expected time. 

CE-certificates. Consider a time to when a CE-certificate at some node v fails. We assume without 



^We emphasize again that arguments of this kind are based on the assumption that the motion of the 
points is oblivious to the choice of priorities. 
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bridge(r, 



left(r^) 




right(T^ 



Figure 6: The view after a swap event between a pair of points p,q with the smallest a;-coordinates stored in 
the left subtree of a node w, whose CE-certificate has to be updated. 



loss of generality that at time to the leftmost vertex p of L(ri;) becomes collinear with the leftmost 
edge qr of 7?."'"(t„), so that bridge(T„) was pq before the event and is pr afterwards, and treat the 
remaining cases symmetrically. See Figure [7] for an illustration. As a result of this event, the edge 
pr replaces pq as bridge(rt,), the edge qr is added to the end of TZ{Ty), and the triangulation T{ti,) 
gains the new triangle Apqr. We need 0{logn) time to update the edge lists of 7^(t„) and T{Ty), 
and to compute the CV-certificate of q (which ceases to be the endpoint of i?(T„)) and add its failure 
time to Q. (Note that the CV-certificate of q is part of the former CE-certificate at u .) 

To recompute the new certificate of bridge(r^), we have to determine the next edge rr~^ of 
7^"'"(r„) that is incident to r from the right. This edge is either stored in one of the lists L(tu,) or 
TZ{tw) at some ancestor w of v, or it belongs to the upper hull of P. See Figure [7] (left). We find 
rr^ by doing a binary search on the lists L(t^) and 'R{tw) for the ancestors w of v, and if necessary 
also on the convex hull of P. 




apex(T„) 



apex(T^) 




apex(r„) 



Figure 7: Failure of the CE-certificate at v (shown at time f^ right after the event). Left: The case where the 
ancestor w that stores rr^ coincides with the ancestor u that has lost q. Right: The case where u is distinct 
from w. 



If pq and qr were part of the upper hull at time t^, we replace them by a single edge pr, in 
O(logn) time. Otherwise, v has some ancestor u such that pq and qr are stored in the edge list 
of L{tu) or TZ{tu). (There is exactly one such ancestor n, which is equal to w unless r is incident 
to bridge(ru); see Figure El In the terminology of Section [3l r^ experiences an envelope event at 
time tQ.) We find u in 0(log n) expected time by searching the edge lists stored at all ancestor 
nodes of v, whose expected number is bounded by O(logn). We then replace pq and qr by pr in 
the edge list of the respective chain L(r„) or TZ{tu), and remove from Q the failure time of the 
CV-certificate of q (within r„). Moreover, we have to retriangulate a suitable sub-pseudo-triangle 
tq of Tu whose boundary, according to Proposition 13.21 has expected complexity O(logre) (see also 
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Figure[5|). To do so, we first determine tq, by locating the edge qv*{q) in the edge list of T{tu), and 
then looking for the lowest (resp., highest) edge above (resp., below) qi'*{q) which is generated by 
a vertex whose priority is smaller than priority (g). We then recursively triangulate tq, as described 
in the proof of Theorem 14.11 All this can be done in 0(log n) expected time. Therefore, we can 
process any CE-certificate failure in 0(log n) expected time. 

Cy -certificates. We finally consider the case when a visibility event, involving some point p within 
the funnel of r^,, for some node v of VT~^{P), causes the failure of the corresponding CV-certificate 
at some time to- Since the failed certificate is associated with an internal vertex of L(t„) or TZ{Ty), 
all the necessary updates are local to the funnel of r.u, and to its triangulation T{tv). We update 
the CV-certificate of p and insert its new failure time into Q, in O(logn) time (the new neighbor 
of ^{p) is easily obtained in O(logn) from the respective edge list). In addition, we may have to 
determine and re-triangulate a suitable sub-pseudo-triangle tq of r^,, whose boundary has expected 
complexity O(logn) (see Proposition 13. 2p . An in the case of a failure of a CE-certificate, this can 
be done in O(log^n) expected time, by searching the edge list of T{Ty). 
We thus obtain our main theorem. 

Theorem 4.2. Let P{t) he a collection of n moving points, as above. We can maintain the 
triangulation T{P{t)) under motion in a kinetic data structure of linear size, which processes an 
expected number o/ 0(?7-^/3s_|_2(n) logn) events, each in 0(log n) expected time, where s is the 
maximum number of times at which any single triple of points of P{t) can become collinear. 

Enforcing locality. As implied by Theorem 14.21 the proposed data structure for maintaining 
T{P) is compact, efficient, and responsive (where the last two properties hold in expectation). 
To make it also local (in expectation), it is sufficient to ensure that at any moment of time the 
expected number of kinetic certificates involving any single point is O(logn). Clearly, each point 
is associated with at most two CT-certificates. Since the expected depth of VT^{P) is O(logn) 
and each pseudo-triangle of VT^{P) defines a single CE-certificate, each point participates in an 
expected number of O(logn) CE-certificates. 

We next slightly modify the definition of CV-certificates, in order to ensure that at any moment 
of time the total expected number of CV-certificates involving any point is also O(logn). Consider 
a fixed moment of time to and a fixed node v in VT^{P(to)), and choose any vertex p on, say, 
the left chain L(Tt,). Currently, p participates in a single certificate that it generates (ensuring the 
validity of i^{p)), and in an arbitrary number of certificates generated by all the vertices q of TZ{Ty) 
satisfying i/(g) = p. We modify our algorithm by keeping (i.e., storing in Q the failure times of) 
only the certificates of p that are generated by the leftmost and the rightmost such vertices q in 
TZ{Ty). If p lies on 7^(r^), we act symmetrically. We apply this modification to every node v and 
every vertex of L(r^,) U TZir^). This modification does not affect the correctness of the kinetic data 
structure because, as can be easily checked, among all the CV-certificates involving p and points q 
with u{q) = p, the first to fail must be the extreme ones that we keep. 

Now, at each node v, every vertex of L(t„),7^(t„) participates in at most three CV-certificates. 
Since the expected depth of VT^{P) is O(logn), the asserted (expected) locality bound follows. 
The kinetic maintenance of this restricted set of CV-certificates resembles that of the original set, 
with the following minor modification. Each time when we process a visibility event caused by the 
failure of some CV-certificate, generated by a vertex p at some node v, we also have to recompute 
the CV-certificates involving the old and the new points ^{p). This can be done in O(logn) time 
using a binary search over L(r.i,) or 7^(r^), which does not affect the time bounds in Theorem 14. 2[ 
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